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ABSTRACT 

Growth  and  development  of  any  nation  is  directly  influenced  by  the  health  of  manufacturing  sector. 
Developments  in  the  manufacturing  industry  serve  as  a reliable  index  that  scales  the  development  of  a society.  This  paper 
summarizes  the  results  of  the  research  work  with  an  aim  to  effectively  schedule  job  shops.  The  primary  objective  of  this 
work  is  to  use  meta-heuristics  as  a tool  for  optimizing  the  scheduling  problem.  In  order  to  present  the  research  to  the  user 
in  the  form  of  a tool  a Graphical  User  Inter  phase  (GUI)  is  designed  and  presented.  This  tool  provides  the  user  with  ease  of 
use  and  flexibility  in  choosing  a particular  schedule.  The  GUI  is  coded  using  MATLAB.  Three  meta-heuristic  approaches 
namely  Particle  Swarm  Optimization  (PSO),  Chaotic  Particle  Swarm  Optimization  (CPSO),  Genetic  Algorithm  (GA)  are 
used  to  optimize  the  objective  function  and  get  an  optimum  schedule.  The  validity  of  the  proposed  methods  is  established 
by  considering  and  analyzing  the  benchmark  problems. 

KEYWORDS:  An  Automated  Tool,  Optimizing  Job  Shop  Scheduling  Using  Metaheuristics,  GUI  is  Coded  Using 
MATLAB 

INTRODUCTION 

Job  shop  scheduling  has  been  proven  to  be  an  intractable  problem  for  analytical  procedures.  Recent  advances  in 
computing  technology,  especially  in  artificial  intelligence,  have  alleviated  this  problem  by  intelligently  restricting  the 
search  space  considered,  thus  opening  the  possibility  of  obtaining  better  results.  Researchers  have  used  various  techniques 
that  were  developed  under  the  general  rubric  of  artificial  intelligence  to  solve  job  shop  scheduling  problems  [1]. 
Scheduling  involves  the  allocation  of  resources  over  a period  of  time  to  perform  a collection  of  tasks  .It  is  a decision- 
making process  that  plays  an  important  role  in  most  manufacturing  and  service  industries  .Scheduling  in  the  context  of 
manufacturing  systems  refers  to  the  determination  of  the  sequence  in  which  jobs  are  to  be  processed  over  the  production 
stages,  followed  by  the  determination  of  the  start-time  and  finish-time  of  processing  of  jobs.  An  effective  schedule  enables 
the  industry  to  utilize  its  resources  effectively  and  attain  the  strategic  objectives  as  reflected  in  its  production  plan. 
The  most  common  manufacturing  system  worldwide  is  the  job  shop. 

Job  shops  are  associated  with  the  production  of  small  volumes/large  variety  products  and  operate  in  make-to-order 
mentions  that  approximately  50  to  75  % of  all  manufactured  components  fall  into  this  category  of  low  volume/high  variety 
and  due  to  the  market  trends  this  percentage  is  likely  to  increase.  Even  though  flexible  manufacturing  systems  are  today’s 
keywords  that  frequently  appear  in  many  research  agendas,  scheduling  of  job  shops  still  receive  ample  attention  from  both 
researchers  and  practitioners  due  to  the  reason  that  job  shop  scheduling  problems  exist  in  many  forms  in  most  of  the 
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advanced  manufacturing  systems.  Besides,  analysis  of  job  shop  scheduling  problems  provides  important  insights  into  the 
solution  of  the  scheduling  problems  encountered  in  more  realistic  and  complicated  systems  [2] . 

Scheduling  has  been  a subject  of  a significant  amount  of  literature  in  the  operations  research  field  since  the  early 
1950s  [2]  [3]. The  main  objective  of  scheduling  is  an  efficient  allocation  of  shared  resources  over  time  to  competing 
activities.  Emphasis  has  been  on  investigating  machine  scheduling  problems  where  jobs  represent  activities  and  machines 
represent  resources.  The  problem  is  not  only  NP  -hard,  but  also  has  a well-earned  reputation  of  being  one  of  the  most 
computationally  difficult  combinatorial  optimization  problems  considered  to  date.  This  intractability  is  one  of  the  reasons 
why  the  problem  has  been  so  widely  studied.  The  problem  was  initially  tackled  by  “exact  methods”  such  as  the  branch  and 
bound  method  (BAB),  which  is  based  on  the  exhaustive  enumeration  of  a restricted  region  of  solutions  containing  exact 
optimal  solutions.  Exact  methods  are  theoretically  important  and  have  been  successfully  applied  to  benchmark  problems, 
but  sometimes  they  are  quite  time  consuming  even  for  moderate-scale  problems.  With  a rapid  progress  in  computer 
technology,  it  has  become  even  more  important  to  find  practically  acceptable  solutions  by  “approximation  methods” 
especially  for  large-scale  problems  within  a limited  amount  of  time  [3].  Stochastic  local  search  methods  are  such 
approximation  methods  for  combinatorial  optimization.  They  provide  robust  approaches  to  obtain  high  quality  solutions  to 
problems  of  realistic  sizes  in  reasonable  amount  of  time.  Some  of  stochastic  local  search  methods  are  proposed  in  analogies 
with  the  processes  in  nature,  such  as  statistical  physics  and  biological  evolution,  and  others  are  proposed  in  the  artificial 
intelligence  contexts.  They  often  work  as  an  iterative  master  process  that  guides  and  modifies  the  operations  of  subordinate 
heuristics;  thus  they  are  also  called  metaheuristics. 

Metaheuristics  have  been  applied  to  wide  variety  of  combinatorial  optimization  problems  with  great  success. 
Particle  swarm  optimization  was  developed  by  Kennedy  and  Eberhart  (1995)  as  a stochastic  optimization  algorithm  based 
on  social  simulation  models.  The  algorithm  employs  a population  of  search  points  that  moves  stochastically  in  the  search 
space.  Concurrently,  the  best  position  ever  attained  by  each  individual,  also  called  its  experience,  is  retained  in  memory. 
This  experience  is  then  communicated  to  a part  or  to  the  whole  population,  biasing  its  movement  towards  the  most 
promising  regions  detected  so  far.  The  communication  scheme  is  determined  by  a fixed  or  adaptive  social  network  that 
plays  a crucial  role  as  to  the  convergence  properties  of  the  algorithm  [4].  The  PSO  algorithm  searches  for  the  best  solution 
over  the  complex  space  through  co-operation  and  competition.  First  of  all,  the  PSO  algorithm  creates  the  initial  particle 
swarm,  namely,  it  initializes  a swarm  of  particle  randomly  in  the  available  solution  space,  making  each  particle  an 
available  solution  of  the  optimization  problem.  Furthermore,  the  target  function  determines  the  fitness  value  through  the 
target  function.  Each  particle  will  move  in  the  space  of  the  solution,  with  its  direction  and  distance  determined  by  speed. 
The  general  particle  will  move  following  the  best  current  particle,  obtaining  the  best  solution  by  searching  generation  by 
generation  [5]  [6].  In  each  generation,  the  particle  will  trace  two  limited  values,  one  of  which  is  the  best  solution,  pbest , 
which  is  found  so  far  by  the  particle  itself.  The  other  is  the  best  solution,  gbest , which  has  been  found  so  far  by  general 
group  swarm  [7].  Genetic  algorithms  belong  to  the  class  of  evolutionary  algorithms.  These  are  algorithms  which  are  based 
on  the  principles  of  natural  evolution,  and  they  can  be  divided  into  four  major  types  of  algorithms:  genetic  algorithms 
(GA),  genetic  programming,  evolution  strategies  and  evolutionary  programming.  All  these  types  of  algorithms  are  based 
on  a population  of  individuals.  Evolutionary  algorithms  have  been  applied  to  many  problems  in  management,  e.g.,  to 
location,  inventory,  production,  scheduling,  distribution  or  timetabling  problems.  The  use  of  evolutionary  algorithms  for 
shop  scheduling  problems  started  around  1980.  Two  of  the  first  applications  to  flow  shop  scheduling  problems  have  been 
given  by  Werner  [8,  9],  and  the  first  application  to  job  shop  scheduling  problems  can  be  found  in  [10].  Genetic  algorithms 
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are  the  most  popular  variant  of  evolutionary  algorithms 

JSP  MATHEMATICAL  MODEL 

Shop  scheduling  problems  belong  to  the  class  of  multi-stage  scheduling  problems,  where  each  job  consists  of  a set 
of  operations.  Among  the  shop  scheduling  problems,  there  are  three  basic  types:  a flow-shop,  a job  shop  and  an  open-shop. 
In  a job  shop  problem,  a specific  technological  route  is  given  for  each  job.  Scheduling  is  the  allocation  of  shared  resources 
over  time  to  competing  activities.  The  n*m  jobshop  scheduling  problem,  designated  by  the  symbols  n/m/G/Cmax  can  be 
described  by  a set  of  n jobs  {ji}  i<j<n  which  is  to  be  processed  on  m machines  {mr}1<r<m  Job  shop  scheduling  problem  shall 
meet  the  following  constrained  conditions 

• Each  machine  can  only  process  one  working  procedure  of  certain  work  piece  in  a period  of  time. 

• Each  working  procedure  will  not  be  interrupted  by  other  working  procedures  during  the  processing. 

• Each  work  piece  shall  experience  the  processing  on  m machines,  and  during  the  processing,  new  work  piece  shall 
not  be  added,  and  the  processing  cannot  be  terminated. 

The  objective  of  optimizing  the  problem  is  to  find  a schedule  that  minimizes  Cmax.  .The  objective  function  is  defined 
as 


Cjr  —min  (c  max) 

{Cjrj\<j<n\\<r<m  is  the  schedule  of  completion  times  for  each  operation  that  satisfies  above  constraints. Cmax  is 
the  time  required  to  complete  all  the  jobs  is  called  the  makespan,  where 

Cmax  = max  1 <j<n;  1 <r<mcjr. 

The  processing  of  job  Jj  on  machine  Mr  is  called  the  operation  Ojr.  Operation  Ojr  requires  the  exclusive  use  of 
Mr  for  an  uninterrupted  duration  pjr , its  processing  time;  the  preemption  is  not  allowed.The  starting  time  and  the 
completion  time  of  an  operation  Ojr  is  denoted  as  s jr  and  cjr 

Respectively.  The  predefined  technological  sequence  of  each  job  can  be  given  collectively  as  a matrix  {T  jk}  in 
which  T jk  - r corresponds  to  the  k-th  operation  Ojr  of  job  Ji  on  machine  Mr. 

PSO  AND  CPSO 

Particle  swarm  optimization  (PSO)  is  a swarm  intelligence  algorithm  that  was  put  forward  through  the  study  on 
the  bird  swarm’s  flying  behaviors.  Similar  to  other  optimization  algorithm  ideology,  one  particle  denotes  one  bird  in  PSO, 
and  each  particle  is  assigned  an  initial  location  and  speed.  During  the  particle  swarm  flying  process,  flying  speed  and 
orientation  will  be  constantly  adjusted,  so  as  to  find  the  optimal  solution  [11].  In  PSO  algorithm,  particle  constantly 
updates  its  speed  and  location  according  to  best  P and  best  g . When  one  particle  finds  one  optimal  local  solution,  other 
particles  will  be  attracted  by  the  optimal  solution  to  gather  around  the  solution  rapidly.  That  will  lead  to  premature 
convergence  and  local  optimization,  which  will  accordingly  influence  PSO’s  search  performance  [12] [13].  Each  particle 
updates  its  position  based  upon  its  own  best  position,  global  best  position  among  particles  and  its  previous  velocity  vector 
according  to  the  following  equations: 
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vfA  =wxv]c+q  xt-  x(j^  —xf)+q  xr2  -g) 


(1) 


_Y* 
*1  ~~  xi 


■Z*Vi 


. A + l 


Where, 

k~ i-i 

P;  : The  velocity  of  Ith  particle  at  (k  + l)th  iteration 
w : Inertia  weight  of  the  particle 

k 

: The  velocity  of  Ith  particle  at  kth  iteration 
c,  c2  : Positive  constants  having  values  between  [0,  2.5] 
rx , r2  : Randomly  generated  numbers  between  [0,  1] 

Pbesti  -The  best  position  of  the  ith  particle  obtained  based  upon  its  own  experience 


(2) 


8 best  : Global  best  position  of  the  particle  in  the  population 
xtk+1 : The  position  of  Ith  particle  at  (k  + l)th  iteration 
x.k  : The  position  of  Ith  particle  at  kth  iteration 
X : Constriction  factor.  It  may  help  insure  convergence. 

Suitable  selection  of  inertia  weight  w provides  good  balance  between  global  and  local  explorations. 

w —w- 

w = wmax ^ — xzter  (3) 

lterm^ 

Where,  wmax  is  the  value  of  inertia  weight  at  the  beginning  of  iterations,  wmin  is  the  value  of  inertia  weight  at  the 
end  of  iterations,  iter  is  the  current  iteration  number  and  itermax  is  the  maximum  number  of  iterations. 

Chaotic  is  a nonlinear  system  that  is  similar  to  the  “Random”  and  has  complex  behaviors.  Since  Chaotic  is 
sensitive  to  the  initial  value,  it  can  easily  jump  out  of  the  local  minimum.  Also,  its  search  speed  is  very  fast.  The  basic 
ideology  for  CPSO  algorithm  is:  In  each  iterative  process,  exert  chaotic  perturbation  on  best  g , and  take  it  as  particle 
updating  position,  so  as  to  prevent  particle  positions  from  converging,  otherwise  it  will  search  locally  around  the  global 
optimal  solution.  A detailed  block  diagram  of  the  proposed  work  is  given  in  the  figure  below. 
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Figure  1:  Block  Diagram  of  the  Proposed  CPSO  Method 
GENETIC  ALGORITHMS 

Genetic  algorithm  (GA)  is  one  of  the  most  widely  used  artificial  intelligent  techniques  for  optimization.  GA  was 
first  developed  by  John  Holland  [14].  GA  is  stochastic  searching  algorithm  based  on  the  mechanisms  of  natural  selection 
and  genetics,  and  is  very  efficient  in  searching  for  global  optimum  solutions.  The  main  idea  of  GA  is  to  mimic  the  natural 
selection  and  the  survival  of  the  fittest  [24].  In  GA,  the  solutions  are  represented  as  chromosomes.  The  chromosomes  are 
evaluated  for  fitness  values  and  they  are  ranked  from  best  to  worst  based  on  fitness  value.  The  process  to  produce  new 
solutions  in  GA  is  mimicking  the  natural  selection  of  living  organisms,  and  this  process  is  accomplished  through  repeated 
applications  of  three  genetic  operators:  selection,  crossover,  and  mutation.  First,  the  better  chromosomes  are  selected  to 
become  parents  to  produce  new  offspring  (new  chromosomes)  [14]  .The  selection  probabilities  are  usually  defined  using 
the  relative  ranking  of  the  fitness  values.  Once  the  parent  chromosomes  are  selected,  the  crossover  operator  combines  the 
chromosomes  of  the  parents  to  produce  new  offspring  (perturbation  of  old  solutions).  Mutation  is  a mechanism  to  inject 
diversity  into  the  population  to  avoid  stagnation.  In  addition  to  the  population  size  and  the  maximum  number  of  iterations, 
several  decisions  on  parameters  must  be  made  for  GA.  Crossover  method  and  crossover  probability  are  the  second  set  of 
decisions  to  be  made.  Finally,  the  mutation  method  and  mutation  probability  must  be  selected  as  they  may  help  to  maintain 
the  diversity  of  the  population  by  injecting  new  elements  into  the  chromosomes.  In  general,  these  three  sets  of  decisions  are 
setempiricallyusingpilotruns.TheflowchartoftheGeneticAlgorithmisinFigure  (2). 
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Figure  2:  The  Flow  Chart  of  the  Genetic  Algorithm  for  the  Proposed  Optimization  Problem 

BENCH  MARK  PROBLEMS 

The  three  well-known  benchmark  problems  with  sizes  of6_6,  10_10  and  20  _ 5 (known  as 

mt06,  mtlO  and  mt20)  formulated  by  Muth  and  Thompson  [15]  are  commonly  used  as  test  beds  to  measure  the 
effectiveness  of  a certain  method.  The  mtlO  and  mt20  problems  are  almost  similar.  They  are  processing  the  same  set  of 
operations  and  technological  sequences  are  similar,  but  in  the  mt20  problem,  the  number  of  machines  available  is  reduced 
to  half  of  that  of  the  mtlO  problem.  For  example,  the  first  operation  of  each  job  in  mtlO  is  exactly  same  as  the  first 
operation  of  each  of  the  first  10  jobs  in  mt20  and  the  second  operation  of  each  job  in  mtlO  is  exactly  same  as  the  first 
operation  of  each  of  the  second  10  jobs  in  mt20.  Taillard  proposed  a set  of  80  JSP  and  120  FSP  benchmark  problems.  They 
cover  various  ranges  of  sizes  and  difficulties.  They  are  randomly  generated  by  a simple  algorithm.  In  this  work  5 
benchmark  problems  are  considered  namely  mt06,  mtlO  and  mt20  formulated  by  Moth  and  Thompson  and  Taillard’ s 15 
jobs  * 15  Machines  and  30  jobs  * 15  Machines  problems. 

MATLAB  GRAPHICAL  USER  INTERFACE  (GUI) 

MATLAB  is  widely  used  in  all  areas  of  applied  mathematics,  in  education  and  research  at  universities,  and  in  the 
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industry.  MATLAB  stands  for  MATrix  LABoratory  and  the  software  is  built  up  around  vectors  and  matrices.  This  makes 
the  software  particularly  useful  for  linear  algebra  but  MATLAB  is  also  a great  tool  for  solving  algebraic  and  differential 
equations  and  for  numerical  integration.  MATLAB  has  powerful  graphic  tools  and  can  produce  nice  pictures  in  both  2D 
and  3D.  It  is  also  a programming  language,  and  is  one  of  the  easiest  programming  languages  for  writing  mathematical 
programs.  MATLAB  also  has  some  tool  boxes  useful  for  signal  processing,  image  processing,  optimization;  etc.MATLAB 
has  many  advantages  compared  to  conventional  computer  languages  for  solving  technical  problems.  MATLAB  is  an 
interactive  system  whose  basic  data  element  is  an  array  that  does  not  require  dimensioning.  The  software  package  has  been 
commercially  available  since  1984  and  is  now  considered  as  a standard  tool  at  most  universities  and  industries  worldwide 

The  GUI  is  coded  using  MATLAB  Version  7.1  and  is  designed  to  enable  the  user  to  have  seamless  analysis  of  the 
data  using  different  optimization  methods.  The  data  required  for  analysis  is  fed  through  an  Excel  sheet  in  predefined 
format.  This  helps  in  standardizing  the  input  methods  and  helps  in  avoiding  user  induced  errors.  A Graphical  User 
Interface  enables  the  user  to  have  seamless  use  and  flexibility  of  operation.  The  implementation  is  carried  out  in  a system 
having  Core  2 Duo  processor  cloaking  at  a speed  of  2 GHz  with  a RAM  of  2GB. 


Figure  3:  Screen  Shot  of  the  Gui  Designed 

From  the  GUI  it  can  be  observed  that  Input  to  the  system  can  be  given  by  pressing  the  Load  Process  Details 
switch.  Once  the  data  is  loaded  the  number  of  machines  and  the  number  of  jobs  involved  in  the  scheduling  are  displayed 
below.  The  user  can  then  choose  the  required  method  for  scheduling.  The  scheduling  results  are  displayed  in  terms  of  make 
span  of  each  proposed  schedule,  a decision  table  which  tabulates  the  job  sequence  and  a Gantt  chart. 

The  functional  icons  present  in  the  GUI  can  be  described  as  below  in  reference  to  the  Figure  3. 

• Functional  icon  used  to  load  the  data,  in  the  form  of  Excel  sheet  which  has  the  machine  sequence  and  timing 
details 

• This  functional  icon  is  used  to  choose  different  optimization  method  for  scheduling 

• The  makespan  of  that  particular  schedule  is  displayed  here 
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• Machine  and  the  job  sequence  is  displayed  here 

• Gantt  chart  of  the  schedule 

RESULTS  AND  CONCLUSIONS 

In  this  work  Genetic  Algorithm  (GA)  function  available  in  the  Matlab  optimization  tool  box  is  used  in  the 
proposed  work.  The  population  size  is  fixed  at  20.  The  elite  count  used  is  fixed  at  10  % of  the  population  which  2.  The 
selection  is  based  on  ranking.  The  cross  over  fraction  is  fixed  at  0.2  and  the  adaptive  feasible  mutation  function  is  used. 
The  migration  of  the  population  is  fixed  as  forward  with  a forward  fraction  of  0.2.  The  maximum  number  of  generations  is 
fixed  at  100.  The  iteration  settings  for  CPSO  include  100  maximum  numbers  of  iterations,  with  acceleration  constant  of  2 
and  2.5  and  maximum  and  minimum  inertia  weights  at  1 and  0.2  respectively.  The  maximum  and  minimum  velocity  of 
particles  is  fixed  at  0.003  and  -0.003  respectively.  The  simulations  are  carried  out  in  a system  having  Core  2 Duo  processor 
cloaking  a speed  of  2 GHz  with  a RAM  of  2GB 


Table  1:  Make  Span  as  Achieved  by  the  Two  Optimization  Methods 


Problem  Size 

Make  Span 
Using  Pso 

Make  Span  Using 
Cpso 

Make  Span  Using 
Ga 

Mt06-  6 jobs  * 6 Machines 

59 

57 

54 

MtlO-  10  jobs  * 10  Machines 

1036 

972 

954 

Mt20-  20  jobs  * 5 Machines 

1328 

1296 

1205 

Tail  lards  15  jobs  * 15  Machines 

1435 

1384 

1317 

Tail  lards  30  jobs  * 15  Machines 

2975 

2765 

2673 

The  Results  of  make  span  achieved  by  different  optimization  methods  of  PSO,  CPSO  and  GA  are  given  in  the 


Table  1 


Figure  4:  Plot  of  Make  span  for  Different  Problem  Sizes  and  Optimization  Methods. 


The  results  are  best  results  achieved  when  each  optimization  method  is  run  for  50  times.  It  can  be  observed  from 
the  above  Table  and  figure  3 the  scheduling  results  produced  by  GA  outperforms  both  PSO  and  CPSO.  When  compared  to 
the  PSO  based  optimization  the  GA  based  optimization  reduces  the  makespan  by  nearly  8 % for  the  Mt06  problem.  There 
is  a similar  reduction  of  8 % for  the  Mt-10  problem  which  comprises  10  machines  and  10  jobs.  There  is  a significant  9.2% 
reduction  in  the  makespan  between  PSO  and  GA  for  Mt  -20  problem.  The  optimization  for  Taillards  problem  also  see  a 
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significant  reduction  in  the  makespan  between  PSO  and  GA. 
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Figure  5:  Gantt  chart  as  Plotted  by  the  GUI  for  Mt  20  Problem 

The  above  figure  gives  the  Gantt  chart  of  the  optimized  schedule  as  suggested  by  genetic  algorithm  for  Mt  20 
problem  in  which  20  jobs  are  processed  across  5 machines. 
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